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CONTRIBUTION TO THE AERODYNAMICS OP ROTAT IN G— WING AIRCRAFT* 

By G. Sis singh 



The conventional calculating methods for rotors of 
rotat ing-wing aircraft with hinged blades are extended and 
refined . 

The chief defect of the investigations up to now was 
the assumption of a more or less arbitrary "mean" drag 
coefficient for a section of the blade. This defect is 
remedied through replacement of the constant coefficient 
by a function of higher order which corresponds to the po- 
lar curve of the employed profile. In that way it is pos- 
sible to extend the theory to include the entire range 
from "autogiro" without power input to the driven "heli- 
copter" with forward-tilted rotor axis. The treatment in- 
cludes the twisted rectangular blade and a nontwisted ta- 
pered blade. 

The investigation is based on Wheatley's report (ref- 
erence 10), which in turn is an extension of the well-known 
studies of Glauert (references 1, 3, 5) and Lock (reference 
2). In the interests of clarity, the principal points of 
view of these analyses will be briefly reviewed. 



SUMMARY AND DISCUSSION 



Proceeding from the air flow and stresses, on a section 
of t?ae blade, the formulas for torque, axial and normal 
thrust of a linearly twisted rectangular blade, and a non- 
twisted tapered blade, are derived. 

The principal advantage of the calculation over ear- 
lier investigations is that it is not necessary to depend 
on a rough estimation of the "mean" drag coefficient when 
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computing a certain rotor as the polar coordinates of the 
blade section in question can he applied direct, independ- 
ent of any assumptions or empirically obtained constants. 

The theory includes all flight stages of rotating- 
wing aircraft from autogiro without power input to heli- 
copter with forward-tilted rotor axis without regular pro- 
peller, and gives results up to coefficients of advance of 
from 0.4 to 0.5 which are in satisfactory agreement with 
wind-tunnel tests. 

t 

The calculation becomes uncertain when, at higher co- 
efficients of advance, the zone of the separated flow with- 
in the rotor area is no longer covered by the stipulated 
assumptions and the effect of the occasional reversed- 
velocity region, for which the substitute functions for 
the air force coefficients on the blade element arc no 
longer applicable, is observed as interference. 

The agreement of the calculation with the experimen- 
tal results will be proved and the effect of twist and ta- 
per on the quality of a rotor discussed in a future article* 



I. NOTATION 



1. Rotor 



P .m 2 , swept-disk area, of rotor. 

S m, rotor radius. 

z - number of blades. 

t m, blade chord = t G (l + p x). 

t Q in, blade chord for r = 0, if the blade form is 

lengthened as far as the axis of rotation. 



p _. taper factor for a tapered blade. 

zt 0 

a - solidity = • 

c~ ' t 0 R 4 p 

r - blade mass constant = — — 
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I mkg s a , mass moment of inertia of rotor "blade 

H 

I = - I ■?, \ r (r-o R) dr 
gj (r) 

eR 

For e = 0, I is identical with the mass mo- 
ment of inertia of the rotor "blade ahout the 
horizontal hinge. 

P( r )Jcg/m, "blade weight per unit of length as f(r). 

Mq. mkg, "blade-weight moment referred to flapping hinge. 

.2. Nondimensional Distances 
Referred to Radius 



Sr. 



sa 



distance of a "blade section at distance 
from the axis of rotation. 



e - distance of flapping hinge from rotor axis. 

B - factor allowing for thrust decrease at "blade ui, 

tip for finite "blade numoer, J 

t 0 (1+0.7 p) 

~B Z 1 — • , i0 

- 1.5 R 

"b - factor allowing for separation flow on return- 

ing "blade at higher coefficients of advance. 



3. Rotor Coefficients 

axial thrust 1 . „„„ 4 . >v „ 

}■ propeller-axis system. 



k nr , - normal thrust ■* 

k„ - vertical thrust "\ . , . 

ST Wind-axis system. 

k g jj - horizontal thrust J 

k^ - torque - 

positive for driven rotor, 

negative when operating as \7indmill; 

forces = coefficient ? U 3 p/2, 
moment = coefficient PR U a p/2, 
p air density. 
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4. Speeds • 

u ra/s, tip of "blade section at distance x. 
U m/s, tip speed at tip .area. 
V m/s, flying speed. 

v^ ra/s, axial flow velocity in the normal- plane; 

v,3_ = V sin a - w, positive if flow is upward. 

w ra/s, rotor induced velocity. 

v m/s, resultant relative inflow at a blade section 
with component; 

U v x in the normal plane, 

U Vy perpendicular to the normal plane. 

oj s** 1 , rotor angular velocity. 

|J. - coefficient of advance. 

"K^ - axial flow 7\<j_ = v^/U. 

5. Force Coefficients 

c„ - lift coefficient, 

c a ~ c a ^-p • 

c a ' - lift gradient = 5.6. 

c - constant lift coefficient for separated flow 

ao 

zones. 

c w - drag coefficient 



w 



( C 0 + C X % + C 2 tt p 3 ) 



c„,c, ,c„ profile constants. 

O 1 d 



factor allowing for 'shifting of blade at higher 
coefficients of advance;' i = 1 for snail coef- 
ficients of advance. 
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coefficient for tangential force 
Cf. = c w cos cp — . Cg_- sin cp. . 

For computing c a 1 and c 0 to c s from the 
profile polar, use an aspect ratio of 
t 0 (1 + 0.7 p) t 
1.5 R 

as "basis. 

6, Angles 

rotor angle of attack (normal plane) to flow. 

blade angle of attack, aerodynamic setting of a 
blade section to normal plane; 

* = * 0 + x *i. 

angle of flow, angle of blade inflow to normal 
plane; tan cp = Vy/ v x . 

operating angle of a blade section at ct^ = cp + 
blade angle of rotation, neutral position back, 
flapping angle referred to normal plane 
P = a 0 - a x cos \|/ - bi sin \|/. 
7. Abbreviations 



b n - n, n 



f (x,\J/) dx d\j/ = / / f (x,\(/) dx d\|/ - 




0 o -p, sin^ 2tt 



o IT 
B 2tt 

f (x,v|/) dx d^f = / / f(x,\|/)' dx d\J/ - 
o V -p. sin \}/ 2tt 



y° ^y* f (x,\i/) dx di 



2 y ^ f(x,\)/) dx dt. 

0 TT 
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II. INTRODUCTION 
1. Plow through Rotor 

If the normal plane of a' rotor faces the air flow at 
an angle a, the axial flow through the rotor is, accord- 
ing to figure 1 : 

v^ = V sin a - w 

where w is the rotor induced velocity which for coeffi- 
cients of advance of around |J. > 0.2 within the swept 
rotor-disk area may "be considered as "being constant. It 
is usually figured at half the amount of that computed "by 
the momentum theory for flow at infinity "behind the rotor. 

According to Kussner's hypothesis (reference 28), the 
comprised air quantity is the quantity of air passing 
through a sphere circumscribing the rotor. On these prem- 
ises tho rotor angle of attack cc is: 

tan a = + - U> 



The assumption of constant axial flow within the normal 
plane is no longer justified for the hovering stage. In 
this case tho actual conditions are "better represented "by 
a distribution of the induced velocity increasing linear- 
ly with the radius. Thereby it can be assumed that the 
mean value w from momentum theory is reached by the par- 
ticular blade elements at distance 0.7 R. (See Isacco, 
reference 19 . ) 

2. Flow on a Section of the Blade 

The equations for the components of the relative flow 
of a blade section at distance x are: 

U v x = U (x + U- sin \|/) in normal plane direction (2) 

U v y = U ^4- 2 || •* M> 3 cos \|/) perpendicular to 

the normal plane (fig. 2) (3) 
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Radial velocities are temporarily disregarded. In 
equation (3) 0 is the flapping angle in relation to the 
normal plane. Its representation is limited to a Fourier 
series of first order of the form: 

B = a 0 - a x cos \j/ - b x sin \}/ (4) 

since the calculation of different examples proved that 
the. terms of higher order had practically no effect on 
the rotor quantities in question. 



3. Air-Force Coefficients at Elade Section 
Expressed "oy Substitute Functions .» 

The significant forces on a section of the blade are 
the lift, the drag, and the tangential force in circum- 
ferential direction. In compliance with modern autogiro 
practice, the center of pressure of the air force is made 
to coincide with the centroid of the blade. As a result, 
the blades arc, while revolving, not subject to period 
distortions conditioned by changing air-force moments. 

The lift coefficient is, as usual, expressed by 

c a ~ c a ' a p (5 } 

The lift gradient c a ' for normal airfoil sections 
may be put at 5.6. Deviating from previous investigations, 
the drag coefficient is replaced by the following function 
of the operating angle 

c„ = 



w = c 0 + c x cx^ + c 2 o^ 2 (6) 



The tangential force coefficient is 

c\ = c w cos cp - c a sxn cp 

In view of the smallness of cp , c^. can be expressed suf- 
ficiently exact with 

c t = c w - c a W 

The sign - in accordance with that for the torque - is 
herewith so chosen that a section of the blade actuated by 
the air forces (autorotation range) corresponds to a neg- 
ative c.|. or , respectively. 
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With cp = a p - d and insertion of equations (5) and 

(6) in equation (7), the substitute function for the tan- 
gential force coefficient follo\vs as 

c t = c 0 + a p (c x + * c a ») + a p 2 (c 3 - c a • ) (8) 

A comparison of the substitute functions with the ac- 
tual variation of. the particular coefficient indicates 
that c w is very well represented within the range of ac- 
tual practice. For c a , on the other hand, the break- 
away of the flow is not taken into account. Consequently, 
the substitute function for c t holds only for operating 

angles up to Op = 15°. 

A survey over the fundamental course of c^. is ob- 
tained from figure 3, where c t with its substitute func- 
tion has been plotted against Op for various blade an- 
gles of attack £. 

The curves disclose the fact known from Schrenk's 
presentation (reference 24) that the a p range of propel- 
ling air forces decreases with increasing blade angle of 
attack d until finally no propulsion is at all possible 

on the particular profile at $ ~ 13°. 

They also confirm, as previously mentioned, the good 
agreement obtaining between the substitute function and 
the actual course up to ct^ i 15°. Then, of course, the 
representation is fundamentally erroneous. 

However, since in normal horizontal flight the cen- 
tral and outer blade elements are st.ruck at small operat- 
ing angles, the representation of the coefficients through 
the chosen substitute functions is admissible. Where, 
particularly at high coefficients of advance, the assump- 
tions are no longer complied with, correction factors are 
empl oyed . * 



The substitution of c a and c u with functions of 

higher order is fairly ohvious. But then complications 
occur in the subsequent integrations according to the 
study. 
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Re ve rsed flow. - ' Sonf ormable to Wheatley (reference 
-10), -the lift in the reversed- velocity, region is written 
in equation (5) as 



which gives for c a 



_ V 7 
v x 



c n = c 



' (- II - ^) (9) 
\ v v / 



a ^ 

For c w we introduce in equation (S): 



0 - .Iz + 



v x 



Herewith 



c w - c 0 



+ c , (- + $) + c 2 (- -2 + (10) 
1 v v v > V v T , 



III. H0V3RING STAGE 



Axial flow and velocitie s_.at a section of the blade . - 

According to elementary jet theory, the mean axial flow 
i s : 



a mean 



'The assumption = - Vy/v x + $ deviating from equation 

(9) for the drag coefficient, affords several simplifica- 
tions in the subsequent calculation. Since the drag in 
reverse flow (a^, = 180°) is higher than in flow from the 
leading to the trailing edge, c v , would have to be re- 
placed \ij a function with other coefficients c 0 to C2 • 
But the coefficients of advance p. < 0.5 in question do 
not justify this action, since the effect of the reversecl- 
velocity region on the forces and moments of the whole r'o— 
tor is quite small. A certain balance is established by 
the fact that the expression - Vy/ v x is, in most cases, 

positive, and then a greater operating .angle - that is, 
higher drag coefficient, in involved. 
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Assuming a linear distribution of the rotor-induced 

velocity over the radius, whereby the above 

" ' ■ * a mean 

is reached at x = 0.7,. it is: 

Zy = x \l x 

with 

d i 1.4 V ' 



Since 



= x for a- = 0, the angle of flow cp 



(fig. 2) is equal to A^ independent of the position of 
the section of the "blade. 



Axial thrust.- A z-blade rot^r with blade chord t 



has an axial thrust of r^. 




/ BR f 

1 i ' O ✓ 

■ y ^ \ * c a" dr 



The factor B allows for the thrust decrease at the blade 
tip, due to the finite blade number. Substitution in the 
conventional manner of the component U v x in the normal 
plane for the resultant inflow velocity v, gives the ax- 
ial thrust coefficient o., 

b y 



r 



sa 



J 



'x 



'a 



dx 



Expressing the blade form and twist by functions 

t = t 0 (1 + p x) 
■b = + x d. 



affords 



B 



• sa 



= a c 



a 



x 2 (A di + € Q + x * 1 ) dx 



(12) 
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The solution is: 



sa 



= o- c a* [| B 3 X dl + | B 3 * 0 - + | B 4 « 4 

+ .p (| B 4 * 0 + | B 4 X dx + | B»« x .) 



(13) 



Torque In accordance with the preceding considera- 
tions, it is: 



-a = u y v x 2 x c w dx + p J v x 2 X 2 c w dx 

o 

- C » / V 3 X ZZ flZ + d Q + X 

a J 1 v x W V 



dx 



B 



- p c. 



2 v 3 



V~ e X 



Iz ( Iz + * . + 



v x ^ v x 



*i )dx] 



Introducing c w according to equation (6) gives, with 



v x = 



f d 
a 



x and v y / v x = A d : 



1 

= y x 3 |c 0 + Cl ( A di + ^ 0 +x * 1 ) 

+ P a (^d 1 + *o + x *i) a }ax 

1 

+ p J X* {c 0 + Cl (* dl +« 0 +X *1> 

+ C 3 (X di +^ 0 H-X * x ) B }dX 

B 

c 1 / x 3 h A (X, + £„+x £ n ) dx 



p c a « /x 4 X di (A dl +^ D +x 



dx 
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The result is: 



k d = [j °o + «1 (| ^ + J* 0 * J *i) 



,1 0 2 1 .2 '<L , 



•0 



+ °a (I V + f *d x *o + I \l x *i 
+ i * 0 a + i * x 8 + §* 0 * x )} 

- «a'U B4 V + S4 V I *x)} 



IV. HORIZONTAL FLIGHT 



Following the study of hovering the true flight stages 
will be treated; that is, those stages in which the rotor 
flow is approximately in direction of the normal plane.. 
They are largely the operating stages of normal horizontal 
flight . 

The limit toward one side is formed "by the autogiro 
without power input (free-wheeling condition with = 0), 

and on the other side, "by the driven helicopter without 
regular propeller hut forward-tilted rotor axis. Between 
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these two limits lie the numerous possibilities of a "heli- 
cogiro," on which.- the power can -"be. -arbitrarily distributed 
over rotor and regular propeller.. 

The study of the twisted rectangular blade and the 
nontwisted tapered blade is carried .out separately. Torque, 
axial and normal thrust are computed. The change-over 
from these forces to the flow-axis system (lift and drag) 
is effected by simple transformation of the coordinates 
(fig.' 4). 

To characterize the particular blade, the coefficients 
of axial and normal thrust and torque are supplemented "by 
the indices v and t for twist and taper, respectively. 

1. Linearly Twisted Rectangular Blade 

The equations for the coefficients of flapping motion 
and axial thrust are taken from Wheatley's report (refer- 
ence 10) which, at a twist of the form 

i> = * 0 + x 

read as follows: 

a ° = I {§ 33 Ni+o.oa m- 3 A d+ i « o ( 3 4 +M . 2 3 2 - I n 4 ) 

+ ( B . + fll . (15) 

a 2JJ, 



B 4 - J u? B 2 

X { A d (b 2 - I U. 2 ) + * 0 (I B 3 +0.106 |^ 3 )+^ 3 4 } (16) 

Tj, 4jj,_B /l 0.0135 u.3 \ 

1 B a + * a 2 0 U * I (17) 
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The "blade constant V for a Q in equation (15) is a 
measure for the ratio of the air-force moments to the mass- 
force moments: 

v t o c a' r4 P 
I 

with * 

I = | J P^j r (r-e R) dr (1.9) 
eS 

For the first cases of e = 0 (flapping hinge in axis of 
rotation) I is identical with the mass moment of inertia 
of a "blade with respect to the rotor axis 

R 

i r 

I =— P / n r d r 

e=o k J (r) 



Here P is the "blade weight per unit length. Assuming 

P = constant (uniform weight distribution across the "blade 

length), the expression for Y further reduces to 

R t 0 

Y = 3 g p c a 1 -p~ 

With the normal values: 

p = 0.125 kg- s 2 m" 4 
g = 9 . 81 m/ s 3 
c a '= 5.6 

it finally gives for this specific case (e = 0, P = con- 
stant ) 



Y = 20.6 



R t, 



(20) 
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Calculation of torque, 

H 2* 



O 0 
• BS « 



0 0 



-/ /'■'(-S(-t-*-*Ml 

(21) 



The lift contribution having been investigated by Wheatley 
(reference 10), the terms with c w muBt be evaluated. They give ; 

1 n 

0 0 
1 2n 

Jv x 2 - x-c w -dx-dy> 

— ,«siny n 
— /tsiny 2* 

—J Jv x *.x-c w -dx-dy>\ (22) 

0 n 

0 0 

1 2n 



+ c 2 [^- + » 0 +x& 1 J}dxdy, 
0 re 

+ c 2 (- 0„ + * di) 2 } iv] • (23) 



— /isiny i 
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a 



(c 0 + c x • # 0 + «2 • #o 2 ) J J v x • & • dx dip 
+ (<•, & t + 2 d 0 c 2 ) .J J Ua: 2 • a; 2 • da; d ip 
-f- c 2 -^ x 2 J J "a; 2- 3? -dxdip 

+ c 2 \ | u„ 2 • .r • da; d »/' 



1 2j 



v x -v y - x-dxdip 



0 0 

1 2jt 



+ 2 ^ • c 2 J J ife • v y ■ x 2 ■ dx dip] (24) 

o o 

After arrangement of the terms, together with the 
contribution of c a taken from Wheatley^s report, the torque 
of the linearly twisted rectangular blade is: 

^ = (c 0 + Cl • & 0 + * V) (1 + ±f -L ^ 

+ j h (c, + 2 & B ■ c,) + ■c 1 + 2& 0 & 1 - c 2 ) [~ + ^ /*j 

+ ^ 2 "j+ W»i- (§ + re ^ 

+ i" h o-x ( \ — J- /* 2 ) + V | y — "J" /* 2 ) — a ° 6 i 

- {A/ S 2 - j- ^) + * (1 + ^ * ^ 
+ T Bi **+ S^) + %(| * 2 -| ^ 2 ) 

+ ^(| £4 + ^ 2i?2 ) + v(| i?4 + ^^ 52 )} • < 25 > 

Normal thrust .- The normal thrust is built up 
from the contributions of the forces at the blade 
in circumferential direction (normal plane) and the 
lift components created by the flapping motion. 

In the normal plane a section of the blade at 
distance x has at blade setting i|/ a normal thrust of 

v2 | t dr ct sin \|/ 

sv 2 !| t dr sin \|/ (c w - c a sin (p ) 
Considering first the share of c w , we find: 



^ 5,1I = 2^r v x 2 ' c V sin v dxdip 

71 '■o 0 

1 2 M 

-j-J J i^ 2 • c OT • sin ip- dxdip 
— nsrnxp n 

■ — n&'a\p 2n 
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These expressions outwardly have a resemblance 
to equations (22) to (24). There sin substitutes 
for 1 x x and accordingly, we have: 

Ami i = ^ [(c 0 + Ci • *o + Ca • #o 2 ) J J °* ' sin V<ixdy) 

+ ( c i • *i + 2 #o #i • c») J J »e 2 • a: • sin y c^x d y 
+ <?» • #i 2 J J "a 2 • a; 2 • sin y dx d y 
+ ' c 2 J ^ Uj, 2 • sin v da; dy 

1 2x 

-f- ( c i + 2 # 0 • c 2 ) J J "a " u v " s i n V ^ 
o o 

12" 

+ 2 # t • c 2 J J Ua; • v v • x • sin y> dx dtp 

00 ... (27 

The share of c a in the normal plane amounts to: 



K 3n II n _ A 



2 JI 



f+J JN> a *.|S-^ + # 0 +**i) sin ydaidy 
0 0 xx 
B 2n 

S V * ' v 1 (b 8 + *° + Si " V dV 
— usinu- « 

"J J , V-(-fj(-g-*.-^ 1 )sin V ,da:d V 

. . . (281 



0 7Z 

a da 



^-[jj v v *- S my>dxd W 
x 

+ J J # 0 • v w • v v • sin ydxdip 



»].... (29) 

The normal thrust of a blade with setting i|f ii 
for flapping angle /3 and lift A, 

- A 0 cos\|r 

according to which; 

B n 

*»nm = — ^2^" J' t 'ic 2 -(-^ + *o+ x&^0- cosy>dxdy> 

0 0 
B 2n 

+J* Ju^^ + ^o + aj^/S-cos y>-dxdy> 
— psiny it 
— /<siny a« 

H-J' J* l 'a 2 '( — ^ — — aj^j/J'COSvdaidvj 



0 it 



^ J JJ y,t • » v • cos ydxdip 
x 

-f- # 0 J j" Us 2 - cos y>dxdy> 



+M1 



Us • a; • cos 



Vda;dv] 



(30) 



H.A.C.A. Technical Memorandum Ho. 921 
Introducing 0 = a 0 - ai cos^ - di sin^, we find: 

*»«m = — yfy" ' H~ 00 v " ' Vy ' 005 v dx d v 

x 

— OjJJua,- v y - cos 2 yidxdip 

x ■ 

x 

-\- d 0 • a 0 j"J Da; 2 • cos tp • dx d y> 
x 

— #o • % J J v x 2 ■ cos 2 ydxdtp 

x 

— #o * *i J J v x • cos v> • sin y > dxdip 

X 

-)- # x • a 0 J J v x 2 • x - cos y>- dxdy> 
x 

— #1 • a l \ \ "a 2 • # * cos 2 y> - dxdip 



— #1 ' bi JJ u * 2 • x • cos v> • sin y • da; dyj 

The ultimate result is the sum of k sni -e- ksnnf 
as coefficient of normal thrust: 



ksnv — O 



(c 0 + Cl . # 0 + C2 • (y /*+y 

+ ( Cl • +2 c 2 A^) (j <« + ify 
+ y c 2 V + 2 c 2 * {1 ,» A, - a, (| - 1 /^j 
+ (2c 2 # 0 + ci) (yi"^d— y «i+y i" 2 «i) 
+ c 2 {y f* V + y ^ 3 a,'- «. *x (y/* a - j^ f 3 ) 

-^•^(y-y" 2 )-^(i"-M' t3 ) 
+v(f* + B*)j 

-i-^^l? 2 -^ 2 )}] (32) 
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2. Wontwisted, Tapered Blad e 

The blade chord is to be a function of the 
radius of the form 

t=«o (l +.p- *);. . . . . . (33) 

Flapping motion .- The coefficients a 0 , a^, 
of the flapping motion are as usual obtained by 
Fourier division of the thrust moment M a . 

Bfi 

M s = ^v i -r-t-c a -e/2dr 



= ?7 a • B? • t a • c a '- e /2 [ JV- * -{^ + 0„) (1 + p x) dx 
o 

— nslnifi 2--r 
-2§vJ-x{^+^(l + P x)dx^ ■ ■ ■ (34) 



0 

2n 

The manner of writing j indicates that the 

TT 

second expression should be allowed for only for 
the region of the returning blade ( \\> = rr-r 2tt ) . 
From equation (34) follows: 

B B 

R*-t 0 -c a '. e l2 = ) V * ■ v "-*- dx + & » J y » 2 • ^ * <** 

B B 
+ /> j* u-c i>„ • a; 2 • da; + p & 0 J IV s • x e • dx 
0 0 
— tisiny —usiny 
— 2 | j Us ■ v v ■ x ■ dx -f- #„ J y x 2 • x ■ dx 
o o 

+ p § v x ■ v v • x* dx + p& 0 § v x 2 • x 2 da; J 

0 0 71 

(35) 

The harmonic analysis of this expression giveB: 

+ sin v |y p ■»„ • 5 s + 0,053 p*& 0 — i- a* ■ B 4 
+ i ^ r\ d ■ B 2 - 1- p* X d + ~^ a, B* 



+ cos y |— y a 0 B 3 — 0,035 ,u 4 a 0 

+ 5 4 + T ^6 1 B 2 

+ ^BV 2 6 1 )}- • • • (36) 
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Since in this Fourier series the constant term 
with the exception of one factor must be equal to a 0 
and the factors of sin vj/ and cos vj/ must equal sero, 
the coefficients of the flapping motion are as follows: 

a ° = i v [t la BS + 0,08 * h + T*° ( Bi + **— ¥ 



/•CO 2 

(37) 

2/i 



B* — ~ B* + p B*(o,& B* — y ft) 

X (^ 2 -^ /« 2 ) + #o (-| * 3 + 0,106 

+ P^o^ + y^ 3 ^)} (38) 
. / 1 , 1 _ , 0,035 ,\ 

** — i 7 1 r ' 39 > 

B* + ±^ + p{0,8B* + ±v*B} 

Axial thrust .- The supplementary axial thrust in 
relation to the rectangular blade with chord t 0 , is: 



0 0 

B 'In 



+ J Jv x *-x(^ + & 0 }dxdy, 
•j sin y n 
■/< sin \f/ 2n 



0 7T 

X 



~ ^ 2 [Xf ' Vx ' v y' x ^ x ^V+#o J^JVt 2, xdxdip 



(40) 



In conjunction with the axial thrust for the 
rectangular blade according to equation (18), it 
therefore gives for the nontwisted, tapered blade 

h H = a c a ' \j h (i? 2 + y + | ^ a, 

+ P #o ( T ^ + y *V~ gy . . (41) 

Torque ,- Here also the additive terms relative to 
the rectangular blade are considered. The share of c a 
in analogy with the consideration of equation (21) is: 

' . x 

Ak dCa =— P ' a 2 ° a [ J j V • a; 2 • <*•«•• V 

x 

+ # 0 j J f x • • a; 2 da: (2y j (42) 
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and the additive decelerating moment from 
c w is, according to equation (84): 



1 2x 

+ (c 1 + 2& 0 e ! ^ J §v x -v v -x 2 -dxdy>] (43) 
o o 

After combining and arranging the terms, 
the final result for the torque of a nontwisted, 
tapered blade is: 

= (co + »o ■ h + V • c «) (t + T ^ — A" *) 

+ * {< a + 4 + » ** (-2- - 1 

- { V ( j- a- - | *) + A, 0. ( j B» + ^ „•) 

+^(i fl4 +T^-^ Ba )} 

+ P [ (Co + <V#o + c s -*. 2 ) (g+ J/* 2 ) + ^(cx+Z^oC,) 

+ V (I - 9^7 ?) - T " a " h + 1 °» ! 

- ca' { V (| - JL + h &0 \± B ' + ± „ 
- /* a 0 6, • 5* + a,» (~ + -g- fi» ^) 



(44) 



Normal thrust .- The share of c w in the normal 
plane is, conformable to equation (27): 

A k,nl ^^Yn [( c o+ c i'^°+ c *'*o a ) v * m x ' sin V dx dtp 

-f - Cj i>» ! • x • sin y • dx dy> 

1 2.T 

+ {cj -f" 2 #o * Cg) § § v x -d v - x-s'myj dx dy>^ (45) 



o o 



and the share of c a in the normal plane, 
according to equation (29) : 
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A k s „ n = — P 2 °^ Ca [ Jj v v * • x • sin v dx dip 

X 

+ &ojj§ v x'V v - X'Smytdxdyij . . . . (40' 

The increment from the lift component 
i6 in accord with equation (31) : 

A k anm = — P *3 W C ° [+ oo J J v x • v v • x • cos y> dx dip 

X 

— J J v x • Vy • x • cos 2 ip dx dip 

x 

— b x v x ■ Vy • x • cos ip sin ip dx dip 

x 

+ #„ a 0 v^-x- cos ip dx dip 
x 

— #o °i v x z - x- cos 2 ipdxdip 

x 

— *o &i JJ «« 2 • « • cos y • sin ip dx dipj (47) 

^7 combining these terms with equation (32) , the 
normal thrust of a tapered blade follows as 



+ (2c a # 0 + c x ) (y/iA,, — i- ai + -i-^ai) 

+ c 2 {|- /* A, 2 + -i- a. 2 - a 0 6 t (-1 - j±- ^ 

- ai (t _ T - a * (t " _ W ^ 

+ ^(l" + i" 3 )} 

+ P [(Co + Cl • 0 O + c 2 V) (y /* + ^i" 4 ) 

+ (2c 1 # 0 + c I )(i-M^-y«i + ^/« I «i) 
+ C2 {A *** + ife"' *i (^ 2 -i" 4 ) 

(fit — l^ 2 )}] (48) 



+ 
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V. CORRECTIONS 

1. Allowance for Finite Flapping-Hinge Distance 
from the Rotor Axis 



So far it had "been assumed that the flapping hinge of 
the "blade coincided with the axis of rotation. But in the 
majority of constructions this hinge is a certain distance 
away from the rotor axis and for such cases the flapping 
motion and, consequently, the aerodynamic characteristics 
of a rotor are erroneously reproduced. 

If e indicates the nondimensional distance of the 
flapping hinge, the thrust moment with respect to this 
hinge i s : 



BR 



eR 



v 2 -x c a t (x - e ) R dr 



B 

r 



B 



r. 



= U S -| t R 2 I / c a t x 2 x dx - e : 
e e 



x c a dx 



(49) 



Since e % 0.1 
integration limit e 
Then 

B 

r 



M g = R' 



XT 



in the practical designs, the lower 
can "be fairly accurately put at 0. 



L J 
o 



c a T x 



B 

x dx - e r 

J 



c a dx 



(50) 



where the first integral expresses the lift moment con- 
formable to the cited considerations with e = 0; the sec- 
ond term comprises the change due to the outwardly shifted 
flapping hinge. 

a). Linearly twisted r ectan gu lar ~bla.de For this mo- 
ment change the insertion of the substitute function for 
c a according to equations (5) and (9) gives: 
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B 



e r v2 ' /z; 



-H sin \|/ 2TT 
- 2e f v x "<g + * Q + z ^) dx 



(51) 



TT 



2ti 

where j again indicates that the second integral is to 

TT 

"be considered only for - tt to 2tt. Neglecting this cor- 
rection, which in fact is unimportant for the final result 
at the practical coefficients of advance, the harmonic 
analysis of the expression for the moment change conforma- 
ble to equation (51) gives: 

AM a e ["l ^ ^ 2 , s /I 



§ * ^ t 0 o a » 



[| *d S 2 + £ 0 (| B 3 + | B) 
+ (i B 4 + \ , S 3*) 

+ sin \J/^|j, d 0 3 s +| n di B 3 

- | a x B 3 + n A d B + | n 2 a x b]- 
+ cos \|/ - ^ n a 0 3 2 



+ | hi 3 3 + J (J. 2 hi B J> | (52) 



This affords as new equations for the Fourier coeffi- 
cients of the flapping moments on a twisted rectangular 
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<h =y [y h B* + 0,080 A » k + -j-*„(b 4 + a 2 52 - y) 
.+y*(^+y^ 3 ) 

-«{l^ 2 + *o(y B3 + |^ 5 ) 

+ ^(T B4 + T^ B2 )}]-7& • • • < 53 > 

— 2j 7i V >< (* 2 "T * 2 ) 

B i --0,5n*B* — eijB 3 — /i< i B) 1 v 7 
+ 0„ (-J 5* + 0,106 ?j + fi X 'B* 
— e B + 2&„B* + y*V £ 3 )} . . . (54) 
4^5 (A + 9i??5 . 3 L.) 

B2 + -2^- e (T B + 1 F) For equa- 

The harmonic analysis of this integral gives: tion (56) 

-J 1 q Tg 2 | a 1 1 1 . t 2 B \ seepage 

+ (y#> + 0,08^) 

+ ♦.(^ + ^^-1^)} 
+ sin v |/i #o B 1 — y <H B 3 + fiX d B + -j p* a x • B 
+ P (y /« *o S 3 + 0,053 /** # 0 — ~ a t 



+ 1 /i A d i? 2 - i- ^ 3 A d + i- ^ a, 5 2 )} 

+ y^B 3 — 0,035 ^ a 0 
+ ±b 1 B*+±-^b 1 B*)]\ (57) 



4 ' 8 

The flapping motion coefficients for a tapered 
blade and a flapping hinge shifted out of the axis of 
rotation (see footnote, p. 26). are then as follows: 

<h=l[ (l-p-e) 0A d B 3 + 0,08 fiH a + J#„(iJ 4 + H*B 3 -\ A* 4 )} 
+ p{~- l d £'-0,094 (Sh + 0„ B* + ~fi' 

-4^ 58 W¥ S3+ ^ a *)}]-:fe • < 58 > 

fli = IB . 

(B*-0,5 f i*B*)(l-pe) + p(o,8B !i — ^J" 2 ^ 3 )— e(^B 3 -^B^ 

x ( 52 -y p 2 ) + *o (y- s 3 + o,io6 ^jj (1 _ p . e) 

+ P (b^ 0 +^B 3 X^ — e(2# 0 B* + 2l d B)j. . (59) 

. s «i-4-^(T + a ffv) + '-T-n} - 

"* (I_ 1 ,.)(i. + ?v)+,(w+i^«)-.gB+«g (eo| 
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"blade and a hinge shifted out of the axis of rotation:* 

h) Kont wi sted, tapered "b la de.- The equation for the 
change in lift moment is: 



B 



A M c 



3. li 



*o c a' P/2 



= e 



dx 



o 



+ d 



.) 



dx 



3 



= e 



v x V dx 



r 

+ *o / 



v x B dx 



4-/ 

o 



P 3 

+ p / " 



,, v y x dx+p 



x dx 



(56) 



Equations (54,55) and (59,60) for defining a x and h 1 

are the result of the simplified assumption that the moment 
of the "blade mass forces 

R 

I = ~ / F/x r(r~e P.) dr 

equals the mass moment of inertia: 

H 



eR 

The exact solution gives two equations with unknown a x 

and 'b 1 , which are obtained 'when the nondimen sional factor 
of sin \j/ in the Fourier series for the thrust moment is 

t>i Um~ I) 

equated to — , and that of cos V is equated to 



O A 

6 *o c * ' H 



Itfote: For equations (57) to (50) inclusive, see pa?e 25. 
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The subsequent calculation of the rotor is effected 
with the normal formulas, wherein the coefficients of the 
new flapping motion, according to equations (53) to (55) 
or (58) to (60) are simply inserted. Strictly speaking, 
allowance should also "be made, when determining the flow 
angle cp at the- blade section, that in the expression for 
T„ , according to equation (3), the term due to the flap- 
ping velocity 

dB , v dB 

x ^ in (x - e) ^ 

changes, though this can be ignored for values of e 0.1 



Change in Profile Drag Coefficient Due to Yaw of 
the Blade at Coefficients of Advance 



In the determination of the resultant flow velocity 
of & blade section during rotation, it had thus far been 
assumed that the flow in the normal plane is always per- 
pendicular to the blade - i.e., the velocity component in 
blade direction was ignored. In reality, the resultant 
velocity in the normal piano consists of the circumferen- 
tial velocity of the particular blade element and the fly- 
ing speod V = \x U (fig. 5): 

v , = x U + |j, U 

ro sultant 

hor.ee the blade section is exposed under a changing angle 
8 while rotating. It is: 

ton 8 = — B_£Os_l (61) 
x + |j, sin \|/ 

Since the polar of a profile, especially the drag, 
changes in yaw with the angle of sideslip 8, an average 
value for this S is formed as function of the coeffi- 
cient of advance. In this manner, a comparison with cor- 
responding r-ind- tunnel tests affords a correction factor 
for the drag coefficient which takes care of the yaw of 
the blade at coefficients of advance. 

Allowance for the rating of the blade elements in the 
various settings is carried out in such a way that the in- 
dividual points of the disk area swept by the rotor are 
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classed according to the dynamic pressure under which the 
"blade element . operate s . 3?or blade setting \|/ , it is: 



1 



6 v x 2 dx 



o 

tan 8. 



'mean _ j 



v/ dX 



H cos \|/ (x + (J. sin \J/) dx 
o 



(62) 



2 

(x + n sin \|/) dx 



When forming the integral - over \J/, it should be ob- 
served that for the particular mean value 6 m the absolute 
amount without consideration of the sign is of interest. 
The integral should therefore express only half a rotor. 
Furthermore, 8 changes sign on the returning blade for 
x - - \i sin \[/. Herewith, we have: 



tan 8 m = 
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Prom the presentation of this function in figure 6, 
it is seen that the mean angle of sideslip occurring at co- 
efficients of advance is not insignificant. For n = 0.4, 
for instance, it • amounts to ahout 20°. 

Allowance for thi-s angle of sideslip is now made in 
such a way that the normal polar is computed with i c w 
instead of with c w = c 0 + c x a, p + c c a p 2 . The correction 

factor i is to he taken from corresponding wind-tunnel 
tests as f (8 m ). • 

In earlier aerodynamic studies of rotors, the effec- 
tive flow velocity of a hlade had "been expressed "by the 
component of the velocity perpendicular to the hlade. - 
Hence the calculation was already made with a reduced dy- 
namic pressure. In the evaluation of tests on airfoils in 
yaw, the coefficient of the drag must accordingly "be re- 
ferred to the flow perpendicular to the hlade. 

To my knowledge there are at the present time no ac- 
ceptable test data of this kind for well-rounded "blades 
such as are employed for rotor "olades. For this reason, 
it is very desirahle t<i make such experiments now, in or- 
der that the effect of the yaw, which at high coefficients 
of advance is certainly of no small influence, may he elu- 
cidated. 

3-. Allowance for Break-away of Flow at Returning 
Blade with Coefficients of Advanco 



In the investigations on the representation of the 
air-force coefficients on the hlade element at the "begin- 
ning of the report, it had already teen pointed out that 
the validity of the preceding considerations hinges upon 
the central and outer parts of a. hlade during rotation 
"being subject to small aerodynamic operating angles. 

However, inasmuch as these assumptions are n'o longer 
always complied with at high coefficients of advance, it 
is necessary to discuss the operating angles Op at the 
different points of .the rotor disk more in detail. Con- 
fined to the normal flight stages, where the flow through 
the rotor is approximately in the normal plane, the disk 
area swept "by the "blades can he fundamentally divided into 
three zones of operating angle (fig. 7). 



30 



N.A.C.A. Technical Memorandum No.. 921 



On the side of the. returning hlade (\|/ = 180° to 360°), 
the hlade elements within the disk' area 

. r = - (j.- E "sin v|/ 

are in the reversed velocity region. 

This is followed "by a zone B of separated flow for 
which the representation of the lift coefficient c a ."by 

c a ' cc-p is no longer justified. This critical region can 

he approximately demarcated from the sound flow with r = 
~"h E sin \|/, ' where h is left to he determined. It is 
the distance x at which, for \|/ = 270°, the value 
ap = 15°. is. reached, according to the previous consider- 
ations.' 

For \[/ = 270° , it is 

v y = A d + x ai 

v x = X - \x 

and consequently, 

P x - p. 



Solving this equation with respect to x, it affords with 
cx-o = 0.2 6 

T, = »J°sl 6^_*)_+_ Ad * (64) 
0.26 - a x - * 

For the rest of the principal share A of the disk 
area the angles cc^ are smaller than 15°, which confirms 

the correctness of the theory of the forces at. the "blade 
element. 



It is advisahle, however, to analyze the working angles 
a,p of the different "blade elements for \J/ = 270°, so as 

to gain some opinion ahout the proper selection of h. 
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The correction for the separated flow, (the shaded area) in figure 
7, is applied in such a way that Ca i» replaced for this range by a 
constant value c a0 3£0.5 c^max* Naturally the share already contained in 
the previous equations for k sa , ksn *nd kd must be subtracted again from 
the assumption c a = c a ' ct p . Quite obviously this kind of correction is 
possible only if b does not exceed the value B, but which will under 
normal conditions be rarely the case with the present day customary 
coefficients of advance. 

The signs for the subsequently computed correction terms are so 
chosen that these expressions become additive to the old coefficients. 
To illustrate: if ksa correct iB the coefficient of axial thrust with 
allowance for the separated flow, and Ak 8 a the correction, it is; 

k sa correct " k «* + Ak sa < 65 ) 

a) Linearly twisted rectangular blade . Axial thrust: 

2/r — bsiny 
A k, av = ^c a , J* § v x • dy> dx 
t —ftslnyj 
2n — bsinv 
- c a ' { J JV (~ + #0 + * # ? ) d y> rfsj] 

2/r — bsin Ijl 



2w 

-*«sinv 

2x —bsiny 8* — bsinv 

— #i § § v x- xdy>dx — J* J v x • v y ■ dip dx^ (66) 
n — tisinyi n— /tsinlp 

With the abbreviated method of writing £ n ■ b n - |T n , the result reads, 

A k sav = - a ■ c a ' Q C a - 1 C m) + ^ <h. (C 3 — fV— C i» a ) 

+ (*»-S-)-A(-r £a - ?i "+^ a ) 
WI^'-i^ + SH] • (67) 

Normal thrust .- Similarly, it is for the normal thrust 

2jt — bsinv 

A k snv = ~ [ — c^j" Jv x 2 ■ sin <p sin y> dy> dx 
7i — /tsiny 
2» — 6 sin ijt 

jt — /<sinv 
2* -—bsiny 

o B 2 -sin yidyidx 
-psiny 
2» — bsiny . 

+ J' # 0 * v x ' »» • sin ydyidx 
n — jisiny 
2x — bsiny 
+ J* J*^ x • v x • »y • a; • sin y> d y> ctej- 
n — sin v 
2n — bsinv 



Ctu'^ £v x -v y -sin y> dy> dx^. (68) 



it — psiny 
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This contains only the shares of the forces in the 
normal plane. The correction from the lift component due 
to flapping angle may be neglected, since this effect 
appears at the decisive blade settings of \|/ = 270° 
primarily as lateral force, while the reaction on the 
normal thrust is quite small. On these assumptions the 
normal thrust correction reads, 

-Wif^-ife^ 2 )] (69) 



Torque ■ 



2n — b sin ^ 



+ 



n — Msiny 

2n — bsiny- 

— C "°S S° x * ' X '~v~ ^ dx \ 

t —ft sin y/ 

2 n — b sin y/ 

n — psiny> 
2n — ftsiny 2« — bsiny/ 

#i §v x • v v • x 2 dtp dx + J Jv v 2 ■ x • dtp dx^ (70) 



x— fi sin y 



71 — psiny 



The solution is: 
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b) Non-twisted tapered blade . 
Axial thrust. Conformably to equation (66) 
the axial thrust correction reads; 

. i — jusin y> 

2«r. — ftsiny 
+ J" J* v x -v v - dy> dx 
n —ft sin y> 
2rc — bsiny 

2?r — bsiny 
.+ J* J'«a! , E'»*a:<iv" ia; |j • • • (- 72 ) 

Which, evaluated, gives; 

+ (| p — J- c *») + ^ % (t» - ev-e/i 2 ) 

(73 

Hormal thrust . Conformably to equation (68) it is: 

2« — i>siny> 

7t — ^siny 
2w — b-siny 
+ J j£) w 2 -sinv>(iv'rf* 
w — p sin v 

2 re — bslny 
Zn — b sin ^ 



£. n — u 91U ^ -j 

+ J J » v 2 • x • sin v V [ • ■ ( 74 ) 



« — /tain v 



Hence 



+ a 0 2 .-Lf,. 2 
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-^(^"-i^l- • • (75) 



Torque , Conformably to equation (70) the torque 
correction reads; 

n — it sin v 
2 3t — 6 sin y/ 
+ J J i>„ 2 • x • d y> dx 
a — u sin y 

2rt — & sin v 
+ — Jv x -v v - x*dy>dx 

n — iusin i/t 
2 it — 6 sin v 



a fir — u sill y -. 

+ j JiV.^d V <*4 (76) 
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Figs. 1,2,3,4,5,6,7 




Circumferential 
direction 



Flight direction 



Normal "~-«^ 



Figure 2.- Velocities and angles at 
a section of the blade. 




Figure 1.- Axial flow through 
normal plane. 
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k sv ■» k 8a cos a - k 8n sin d 

k B h = k 8 a g i n <*■ + k sn cos a 
Figure 4.- Air forces and rotor 
setting. 
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Direction of 
rotation 



a. Actual variation 

b, Variation of substitute function 
Figure 3.- Tangential force coeffi- 
cient as f(c6p) for various 

blade angles of attack $ . 
Airfoil section: Gbttingen 367 
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Figure 5.- Blade section in yaw. 





0,¥ 0,6 US ft fff 

Figure 6.- Mean angle of sideslip 

plotted against coef- 
ficient of advance. 



Direction of rotation of rotors 



A, Region of sound flow . 

B, Region of separated flow 

( idealized) 

C, Reversed- velocity region 
Figure 7.- Operating angle of blade 

sections in rotor-swept- 

disk area. 
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